arXiv:astro-ph/0512196v3 20 Jul 2006 


On the “canonical behaviour” of the X-ray afterglows 
of the Gamma Ray Bursts observed with Swift’s XRT 

Shlomo Dado^, Arnon Dar^ and A. De Riijula^ 

ABSTRACT 

The “canonical behavionr” of the early X-ray afterglows of long-duration 
Gamma-Ray Bursts (GRBs) —observed by the X-Ray Telescope of the SWIFT 
satellite— is precisely the one predicted by the Gannonball model of GRBs. 

Subject headings: gamma ray burst 


1. Introduction 

Within a year of its launch on 20 November, 2004, SWIFT accomplished most of its 
goals (Gehrels et ah 2004): it saw and localized short-duration GRBs and discovered their 
X-ray afterglow (e.g. Gehrels et ah 2005; Govino et ah 2005; Better et ah 2005), heralding 
the discovery of their host galaxies in ground-based follow-up observations (e.g. Bloom et 
ah 2005a, 2005b; Hjorth et ah 2005; Antonelli et ah 2005). SWIFT has detected GRBs at 
large redshifts (Jakobsson et ah 2005), with a record 2 ; = 6.29 (Haislip et ah 2005; Gusumano 
et ah 2005; Tagliaferri et ah 2005a). Its X-ray telescope (XRT) measured well the X-ray 
afterglows (AGs) of long-duration GRBs in the first hours after burst (e.g. Ghincarini 2005). 
The latter observations have been claimed to provide two major surprises: 

(a) The 0.2-10 keV light curves of the X-ray AGs exhibit a “canonical behaviour”, to wit: 
(i) an initial very steep decay (f““ with 3 < a < 5), followed by (ii) a shallow decay (0.2 < 
a < 0.8), which finally evolves into (iii) a steeper decay (1 < a < 1.5). These power-law 
segments are separated by the corresponding “ankle” and “break” at times 300 s < fankie < 
500 s and 10^ s < fbreak < 10^ s (Ghincarini et al. 2005; Nousek et al. 2005; Gusumano et 
ah 2005; Hill et al. 2005; Vaughan et al. 2005; Tagliaferri et al. 2005b; Barthelmy et al. 2005). 
This is illustrated in Figs. la,b for GRB 050315 and GRB 050319. 
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(b) While most of the early AG light curves decline smoothly, a substantial fraction has large 
X-ray flares on short time scales (see, e.g. Burrows et ah 2005a, 2005b; Nousek et al. 2005). 

It has also been claimed (e.g. Chincarini et al. 2005; Nousek et al. 2005; Zhang 2005) 
that the “canonical behaviour” and the flares cannot easily be explained by current models. 
This may be true of the popular hreball (FB) models (for a general review, see e.g. Zhang & 
Meszaros 2004). It is not true of the “Cannonball” (CB) model (Dar & De Riijula 2000, 2004 
and references therein). In the CB model the observed canonical behaviour and the flares 
are predictions (Dado, Dar & De Riijula 2002, hereafter DDD2002). This is discussed in 
this letter for two recent representative SWIFT-XRT observations: CRB 050315 (Vaughan 
et al. 2005) and CRB 050319 (Cusumano et al. 2005). 


2. The Cannonball model of GRBs 

In the CB model long-duration GRBs and their AGs are produced by bipolar jets of 
CBs ejected in ordinary core-collapse supernova explosions. An accretion disk or torus 
is hypothesized to be produced around the newly formed compact object, either by stellar 
material originally close to the surface of the imploding core and left behind by the explosion¬ 
generating outgoing shock, or by more distant stellar matter falling back after its passage 
(De Riijula 1987). As observed in microquasars (e.g. Rodriguez & Mirabel 1999), each 
time part of the accretion disk falls abruptly onto the compact object, a pair of CBs made of 
ordinary matter are emitted with large bulk-motion Lorentz factors 7 , in the polar directions 
wherefrom matter has already fallen back owing to the lack of rotational support. 

The 7 -rays of a single pulse in a CRB are produced as a CB coasts through the SN 
“Glory” —the SN light scattered by the SN and pre-SN ejecta. The electrons enclosed in the 
CB Compton up-scatter to CRB energies the photons of the Glory (Dar & De Riijula 2004). 
As the CBs escape the SN environment, the ambient-light distribution becomes increasingly 
thin and radial. No longer can the Glory’s photons be efficiently up-scattered. The 7 -ray 
emission is taken over by thermal bremsstrahlung (TB) and line emission (LE) from the 
rapidly expanding CBs. Their rapid expansion stops (DDD2002) within a few minutes of 
the ejection (of near-axis observer’s time) by their interaction with the interstellar medium 
(ISM). The fast-declining TB and LE is subsequently taken over by synchrotron radiation 
from swept-in ISM electrons spiraling in the CBs’ enclosed magnetic held (DDD2002). The 
CB radiation is a sum of these mechanisms which produce a continuous X-ray AG light 
curve as the dominant CB radiation mechanism evolves according to: Inverse Compton 
Scattering -H’ Bremsstrahlung -|- Line Emission ^ Synchrotron Radiation. 
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In the CB model, the typical parameters of CBs are not predicted but can be extracted 
from the analysis of optical and radio AGs (DDD2002, Dado, Dar & De Riijula 2003a). Their 
values can then be used to predict the properties of the 7 -rays of GRBs (Dar & De Riijula 
2004). They also agree with the characteristic “canonical” timing and flux of the above 
sequence, and the corresponding power-law and spectral changes (DDD2002). The surprise 
(a) of the Introduction was already visible in the X-ray data of GRBs 970508, 970828, 990510 
and 010222 , and compatible with less-precise data on six other GRBs, including 980425. As 
for surprise (b), late-time bumps and flares have been detected before in the X-ray AG of 
GRB 970508, in the radio AG of GRB 030329 and in the optical AG of e.g. GRBs 000301c and 
030329. Their GB-model interpretation is straightforward (DDD2002, Dado et ah 2004b). 


3. Radiation from expanding CBs 


During the GB expansion phase, the emission from a GB is dominated by TB and LE. 
In the rest frame of a GB, interactions of the ISM particles and photons with the GB’s 
constituency produce a quasi-thermal electron population with a power-law tail (oc E~p) of 
Fermi-accelerated (p~2.2) and “Bethe-Block” (p~2.0) knocked-on GB electrons (Dar & De 
Riijula 2001). Thus, a GB emits a TB spectrum with a power-law tail, Doppler-boosted by 
the GB’s motion. The spectral shape (Dar & De Riijula 2004) is: 
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where q;~1,/5 = (p-|-2)/2«i2.1,6isa dimensionless parameter and Tefj = 5T/{l + z), with 
T the GB’s plasma temperature and 5 = 1/(7 (1 — /5 cos 6 ')] 27/(1 -|- 7 ^ 0 ^) the Doppler 

factor of a GB with a bulk Lorentz factor 7 ^ 1 , observed from an angle 0 <C 1 relative to 
its motion. The indexes a and (3 may deviate from the central predictions, as the radiation 
becomes dominated by LE. Also the power-law index of the radiating electrons (after cooling) 
may be larger than p -|- 1 ~ 3.2 (Dar & De Riijula 2004). 


The observed energy flux from a GB at a luminosity distance Dl is: 

^ ^ 3A(r)iv7- 

dt ~ ' ’ 

where Ni, is the GB’s baryon number, R its radius and A(T) its “cooling function”. If the 
loss-rate of the GB’s internal energy is mainly adiabatic, Tocl/i?^. At a GB’s transparency 
time, r, T~ 10^-10® K (DDD2002), and A(T) oscillates and depends on composition in this 
T-range. A rough description of the results of Sutherland & Dopita (1993) is: A(T) ~T“, 
with a ~ 2 for zero metallicity, and a ~ 0 for high metallicity. During the short TB-I-LE 
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phase, 6 stays put and R increases approximately linearly with time. The observer time t is 
related to the CB’s rest-frame time t' throngh dt = (1-1- 2 :) dt'/5. Thus, dF/dt oc 
and the expected powers are~t“^ to~t“'^ (O’Brein et ah 2006 report cases, but they 
were not “canght” early enongh). Here we adopt an intermediate metallicity, a~l, for which 
dF/dt(x(t + t)~^. The CB-model predictions for the identity of the LE-phase lines and the 
time-evolution of their energy are qnite remarkable (Dado, Dar & De Riijula 2003b) —and 
perhaps worth testing. 


4. Synchrotron AG from decelerating CBs 


A CB is assnmed to contain a tangled magnetic held in eqnipartition with the ISM 
protons that enter it. As it ploughs through the ionized ISM, it gathers and scatters its 
constituent protons. The re-emitted protons exert an inward pressnre on the CB, countering 
its expansion. In the approximation of isotropic and complete re-emission in the CB’s rest 
frame and a constant ISM density n, one hnds that within minutes of observer’s time t, a 
CB reaches an asymptotic radius i?( 7 o), with 70 its initial 7 . Subseqnently, ■jit) obeys: 
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with K = 0. If the re-emitted ISM particles are a small fraction of the intercepted ones k = 1 
in Eq. (3). In both cases 7 and 5 change little as long as t < fbreak ~ [l-|-(3 —k) d^ 7 o)] R. 


In the CB model, the ISM electrons that a CB gathers are Fermi-accelerated in the 
CB’s enclosed magnetic maze and cooled by synchrotron radiation to a broken power-law 
distribntion with an injection “bend” at the energy Eb = nieC^'^it) at which they enter 
the CB. Their emitted synchrotron radiation has a broken power-law form with a bend 
freqnency, Vb ~ (1.87 x 10^ Hz) [yit)]^ S{t) [n/(10“^ cm“^)]^/^/(l -|- z), corresponding to Eb- 
In the observer frame, before absorption corrections (DDD2002): 

E^ = u {dn^/dv) oc nR^ [ 7 (^)]^“~^ ['^(^)]^^° > ( 4 ) 


where a ~ 1.1 for u ^ z/^, as in the X-ray domain. The initial slow decline of 7 (t) and 
S{t) resnlts in the shallow decay of the early X-ray synchrotron AG, which is smooth if the 
intercepted ISM density and the extinction along the line of sight are constant. The snm of 
TB and synchrotron emissions produces the canonical X-ray light curve with an early fast 
TB decay, overtaken at the “ankle” by an initially much flatter synchrotron emission, which 
becomes steeper aronnd the “break”, as demonstrated in Figs. la,b. 
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5. Comparison with recent representative observations 

The 0.2-10 keV X-ray AG light curves of GRB 050315 (Vaughan et al. 2005) and GRB 
050319 (Gusumano et ah 2005) —which exhibit the “canonical behaviour”— are compared 
in Figs. la,b, with the sum of the TB-I-LE and synchrotron radiations, Eqs. (2) and (4); 
with use of 'y{t) as in Eq. (3) with k = 1. The two normalizations were best-£t, along with 
the GB model parameters, 70 = 1162, 9 = 0.83 mrad, r = 110 s and to = 0-14 days for GRB 
050315 at z = 1.949, and 70 = 710, 6 = 0.39 mrad, r = 20 s and to = 0.17 days for GRB 
050319 at z = 3.24, all within their usual range (DDD2002). The parameters n, Ni, and 
R occur only in the combination to of Eq- (3) and in the fit normalization of Eq. (4). The 
remaining parameter -the normalization of A(T) in Eq. (2)- occurs only in the corresponding 
flux normalization. Similarly good fits are obtained with 7 (t) as in Eq. (3) with k = 0. More 
accnrate late-time data are reqnired to tell apart the k = 0,1 deceleration laws. 


6. Bumps and flares in the AG 

A fraction of early X-ray AGs have flares on short time scales (see, e.g. Burrows et 
al. 2005a, 2005b). In the GB model, snch early time X-ray flares may be part of a long GRB 
due to late time accretion episodes of the fail-back material on the compact central object, 
and/or if a significant precession of the jet axis occurs prior to these episodes, an X-ray 
flare may be a GRB pnlse viewed more off-axis, with the corresponding time-dilation and 
energy-softening. Another possible sonrce of flares is density jnmps along the GB trajectory 
(e.g. DDD2002; Dado et al. 2004b). In the GB model, the AG is a direct and quasi-local 
tracer of the density of the ISM through which a GB travels; see Eq. (4). Such density jumps 
are produced in the circnmburst environment by the SN and pre-SN ejecta, by the “winds” 
inside the superbubbles (SBs) where most core-collapse SNe take place and, in particular, at 
the complex SB boundaries created by stellar winds and previous SNe within the SB. This 
is demonstrated in Fig. 2a (Dado et al. 2004b) for a density profile shown in Fig. 2b. Fitting 
bnmps to density profiles is not over-informative, unless they are seen at various frequencies, 
as the GB model predicts the energy-dependence of the bnmp widths (Dado et al. 2003a). 
The oscillations of A(T) in the relevant T-range may also indnce X-ray flares. 


7. Discussion 

The early-time X-ray afterglows of GRBs measured with SWIFT-XRT provide an ex¬ 
cellent test for two contenders: the FB and GB models. We have shown that the GB model 
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passes this test with flying colours: the observed “canonical” behaviour was correctly pre¬ 
dicted (DDD2002) long before the launch of SWIFT. To the best of our knowledge, no 
satisfactory FB model explanation of this behaviour has been found, though certain cor¬ 
relations have been proposed. Let the early X-ray flux be described as Fh(t) oc 
The FB-model prediction is ai— (5x = ‘2, as recently discussed by Kumar et al. (2006). The 
CB-model prediction is {3x — 1.1 at all times, so that ai—j3x ~ ai — 1.1. The models are 
compared with the data of O’Brein et al. (2006) in Fig. (3). While the prediction of the CB 
model agrees with the data, the prediction of the FB model is in clear contradiction with 
the observations. 


A celebrated prediction of “conical fireball” models is a break in the AGs when the 
beaming angle of radiation from a decelerating cone increases beyond the opening angle of 
the jet, 7 (^)“^ 9j, and the observer begins to see the full cone (e.g. Rhoads 1999; Sari et 

al. 1999). If the observer happens to be nearly on the cone’s axis, the break time (Sari et 
al. 1999) is: 


^break ~ 2.23 (1 -|- z) 
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with 7]^ the conversion efficiency of the ejecta’s energy into 7 rays. Frail et al. (2001) suggested 
that the 7 -ray energy of conical GRBs is a 6 *j-independent standard candle. Gonsequently 
Eiso ^ the Frail RelationL From 16 GRBs of known ; 2 , Bloom et al. (2003) found 

Ey zz 1.3 X 10^^ ergs. Insert this and 6‘j = 4E^/E1^° into Eq. (5), to obtain: 


fbreak ~ 4.33 (1 z) 
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All published attempts to use Eq. ( 6 ) to predict tbreak —before it was measured— from the 
observed failed^. In fact, in many GRBs the AG is not achromatic as expected in the 
jetted FB model, and the break in the X-ray AG is not matched by a similar break in the 
optical AG (e.g. Panaitescu et al. 2006). This further questions the FB-model interpretation 
of the AG break (Rhoads 1997, 1999; Sari, Piran and Halpern 1999) and the proclaimed 
success of the “Frail Relation”. These suggest that the success of the Frail relation is an 
artefact resulting from an a-posteriori adjustment of free parameters. Moreover, E!^° for all 
XRFs with known 2 ; is much smaller than the “standard-candle” value of Frail et al. 2001, 


^The fraction of the sky illuminated by a conical GRB is /& = (1 — cos 0j)/2 w 6l|/4 and not (1 — cos 9j) w 
d|/2, used by Frail et al. (2001) and many other authors. Bipolar GRBs do light a fraction 6*|/2. 

^For instance, Rhoads et al. 2003 predicted tbreak > 10-8 days for GRB 030226, while Greiner et al. 2003, 
shortly after, observed tbreak ~ 0.8 day. 
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implying that XRFs and GRBs cannot be the same phenomenon viewed from different angles, 
contrary to indications (e.g. Dado et al. 2004a). Many other more successful relations, such 
as the one between the equivalent isotropic energy and the “peak” 7 -ray energy of GRB 
pulses (the “Amati Relation”) are predictions of the GB model (Dar & De Riijula 2004; 
Dado & Dar 2005) but not of the FB model. 

In the GB model, AG flares follow either from GB encounters with density inhomo¬ 
geneities (DDD2002; Dado et al. 2004b) or from late accretion episodes on the compact 
central object. In the FB model, late-time flares result from late central activity (e.g. Gra- 
not, Nakar and Piran 2003). Although such an activity can neither be predicted nor ruled 
out, it is not clear why the ensuing ejecta do not also produce 7 -ray pulses and why the 
duration and magnitude of the AG flares scale roughly with the time and magnitude of the 
declining AG. 


We thank N. Soker for suggesting that LE dominates the GB radiation during the 
TB-I-LE phase. S. Vaughan and G. Gusumano have kindly sent us their tabulated XRT 
measurements of the X-ray AGs of GRB 050315 and GRB 050319, respectively, prior to 
publication. This research was supported in part by the Asher Fund for Space Research at 
the Technion. 
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Fig. 1.— (a) Comparison between the X-ray afterglow (0.2-10 keV) of GRB 050315 measured 
with XRT on board SWIFT (Vaughan et ah 2005) and the CB model £t. (b) The same 
comparison for GRB 050319 (Cusumano et ah 2005). 
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Fig. 2.— (a) Comparison between the R-band AG of GRB 030329, shown as “residua” 
AR of the data relative to a broken power law of index —a jumping from ~ 1.1 to ~ 2 at 
f ~ 5 days (Lipkin et ah 2004), and the residua, relative to the same broken power law, 
calculated from the CB model (red line) for the input density prohle shown in Fig. 2b. The 
K feature is a prediction (Dado et ah 2004b). (b) The density prohle (relative to a smooth 
ISM density —a constant plus a “wind” contribution decreasing as 1/r^). The density is 
n = 'RijTij {Tj/rY 0(r — r^), with 0 Heaviside’s function. 
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Fig. 3.— The CB- and standard (FB-) model predictions for the relation between the indices 
of a simple description of the early X-ray AG flux: F^yit) oc Data analysis from 

O’Brein et al. (2006). 















































